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We show that the non-perturbative vacuum structure associated with neutrino mixing leads to 
a non-zero contribution to the value of the cosmological constant. Such a contribution comes from 
the specific nature of the mixing phenomenon. Its origin is completely different from the one of the 
ordinary contribution of a massive spinor field. We estimate this neutrino mixing contribution by 
using the natural cut-off appearing in the quantum field theory formalism for neutrino mixing and 
"^j ■ oscillation. 
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• By resorting to the recent discovery of the unitary inequivalence between the mass and the flavor vacua for neutrino 
fields in quantum field theory (QFT) [1-8], we show that the non-perturbative vacuum structure associated with 
t-H . neutrino mixing may lead to a non-zero contribution to the value of the cosmological constant. 

The contribution we find comes from the specific nature of the field mixing and is therefore of different origin with 
respect to the ordinary well known vacuum energy contribution of a massive spinor field. 

The nature of the cosmological constant, say A, is one of the most intriguing issues in modern theoretical physics 
and cosmology. Data coming from observations indicate that not only A is different from zero, but it also dominates 
■ the universe dynamics driving an accelerated expansion (see for example [9,10]). 

\ In the classical framework, A can be considered as an intrinsic (i.e. not induced by matter) curvature of space-time 
or a sort of shift in the matter Lagrangian. In the latter case A can be considered as an unclustered, non interacting 
| component of the cosmic fluid with a constant energy density p and the equation of state p — —p. These properties 
^j*. and the fact that the presence of a cosmological constant fluid has to be compatible with the structure formation, 
jjj^ 1 allow to set the upper bound A < 10~ 56 cto -2 [11]. 

In the quantum framework, the standard approach is to consider the cosmological constant as a gravitational effect 
of vacuum energy [10,12-14]. A common problem of all these approaches is that they do not provide a value of A in 
the bound given above. This is known as the cosmological constant problem [10] and it has been tried to address it in 
different ways (see for example [15]). 

In this paper, we show that the vacuum energy induced by the neutrino mixing may contribute to the value of 
cosmological constant in a fundamentally different way from the usual zero-point energy contribution, as already 
stated above. 

Indeed, it has been realized [1,5] that the mixing of massive neutrino fields is a highly non-trivial transformation 
in QFT. The vacuum for neutrinos with definite mass is not invariant under the mixing transformation and in the 
infinite volume limit it is unitarily inequivalent to the vacuum for the neutrino fields with definite flavor number. This 
affects the oscillation formula which turns out to be different from the usual Pontecorvo formula [16] and a number 
of consequences have been discussed [5]. 

The existence of the two inequivalent vacua for the flavor and the mass eigenstate neutrino fields, respectively, is 
crucial in order to obtain a non-zero contribution to the cosmological constant as we show below. 

In Section II, we shortly summarize the QFT formalism for neutrino mixing. In Section III, we show that the 
neutrino contribution to the value of the cosmological constant is non-zero and then we estimate its value by using 
the natural scale of neutrino mixing as cut-off. The result turns out to be compatible with the above mentioned 
upper bound on A. Section IV is devoted to the conclusions. The appendix is devoted to the tetradic formalism in 
the Friedmann Robertson Walker (FRW) space-time. 
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II. NEUTRINO MIXING IN QUANTUM FIELD THEORY 



The main features of the QFT formalism for the neutrino mixing are summarized below. For simplicity we restict 
ourselves to the two flavor case. Extension to three flavors can be found in refs. [8]. The following Lagrangian density 
describes the Dirac neutrino fields with a mixed mass term: 



C(x) = * f (x)(i fl-M)V f (x) . 



(1) 



by 



The relation between Dirac fields &f(x), eigenstates of flavor, and Dirac fields \? m (x), eigenstates of mass, is given 



U is the mixing matrix and = {u\, 1/2)- The mixing matrix is 



U 



cos 9 sin 9 
- sin 9 cos 9 



(2) 



(3) 



being 9 the mixing angle. Using Eq.(3), we diagonalize the quadratic form of Eq.(l), which then reduces to the 
Lagrangian for the Dirac fields ^ m (x), with masses mj, i = 1, 2 : 



where = e?mg(mi, to 2 ). The mixing transformation (2) can be written as [1] 

v a (x) = Gg 1 ^) Vi(x) G e (t), 
where (a, i) = (e, 1), 2), and the generator Gg(t) is given by 



G e (t)=exp 9 J d 3 x (v\{x)v2(x) - v\{x)vi(x)} 



(4) 
(5) 

(6) 



The free fields Vi (i=l,2) are given, in the usual way, in terms of creation and annihilation operators (we use t = xq): 

d 3 k 



(2n)i 



1,2 



(7) 



«w feii t,,r an( j _ ^7^2 _|_ m 2^ rp^g mass eigenstate vacuum is denoted by 



with ul t (t) = e-^ul-, vl-(t) = e-^ V] 
|0) m : a k i \0) m — i\tym — 0- The anticommutation relations, the wave function orthonormality and completeness 
relations are the usual ones (cf. Ref. [1]). 
The flavor fields are obtained from Eq. (5): 



d 3 k 

(2n)i 



«k,i(*)«k lff (*) + «lk,i(*)#l *(*) e *' X ' ^ l ) = ^ 2 )- 



(8) 



The flavor annihilation operators are defined as a k a {t) = G g 1 (t)d^ i Ge(t) and /?l ko .(£) = G 1 (t)(3_ k fig^). They 
annihilate the flavor vacuum |0(t)) / given by 



W))f = G e \t) |0) T 



(9) 



The crucial remark is that, in the infinite volume limit, the vacuum |0(t))/ for the flavor fields and the vacuum 
|0) m for the fields with definite masses are unitarily inequivalent vacua [1]. For a rigorous general proof of such 
inequivalence for any number of generations see Ref. [2]. 

In Ref. [6] it was shown that the correct vacuum to be used in the calculation of oscillation formulas is the flavor 
vacuum |0)j. Indeed, using the mass vacuum |0) m leads to violation of probability conservation [6]. We will thus use 
|0)y in our computations in the following. 

For further reference, it is useful to list the explicit expressions for the flavor annihilation/creation operators. In 
the reference frame k = (0, 0, |k|) they are given by [1]: 
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a r k Jt) = cos0 a k>1 + sin0 ([£(t) a k , 2 + e r F k (i) /3^ 2 ) , (10) 

«£,„(*) = oosfl al i2 - sin0 (l/ k (t) a k>1 - e r T/ k (t) /T^J , (11) 

^ kie (t) = cosfl/3L kil +sin(9(c^(t) / 9L ki2 -e'-V1 c (t)a^ 2 ) , (12) 

^k, M (t)=<^«/^k,2-sin»(^k(t)/3!l k ,i+c r V r k (t)ai;; i ) . (13) 
J7 k and V k are Bogoliubov coefficients given by: 

Vk(i) = |Vk| e i{uk - 3+u *- l)t , y k (t) = |l/ k | e*^*."-"*.!)*, (14) 

l^kl ^ f * f ^1±^\ " U + ^ ) s cos (^), (is) 

it/ i /^M + m A 5 +m 2 \ 1 ( |k| |k| \ . k 

\U k \ 2 + |Vk| 2 = 1. (17) 

The function |Vic| is related to the condensate content of the flavor vacuum [1]: it is indeed zero for mi = m 2 , it 
has a maximum at |k| = y /m 1 m 2 and, for |k| » ^/m 1 m 2 , it goes like | Vk | 2 ~ (m 2 — m 1 ) 2 /(4|k| 2 ). 

III. NEUTRINO MIXING CONTRIBUTION TO THE COSMOLOGICAL CONSTANT 

The connection between the vacuum energy density (p vac } and the cosmological constant A is provided by the well 
known relation 

</W> = (18) 

where G is the gravitational constant. 

In order to calculate {p vac ) we have to consider tetrads and spinorial connection. The symmetries of the cosmological 
metric make this task easier. 

Tetrads are defined as a local inertial coordinate system at every space-time point defined by g^ u — e^e^rjab, where 
g^v it the curved space-time metric and r] a b is the Minkowski metric. In the tetrads framework the parallel transport 
in a torsion free space-time is defined by the affinc spin connection one form lu£ = oo^dx 11 , which satisfies the Cartan 
structure equations 

de a + Lot A e b = 0, (19) 

where e a is the tetrad one-form defined by e a = e^dx v . The energy-momentum tensor is obtained by varying the 
action with respect to the metric g^ u : 



2 SS 
'Sg^(x)' 



= (20) 



where the action is 



S = j ^L(x)d A x. (21) 
In the present case, the energy momentum tensor is given by 

?A X ) = \(^m{x)l^D^ m {x)) , (22) 
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where D v is the covariant derivative: 

D v = d v + T v , Y v = ^ujt[ la , lb \, ltl (x) - 7 c e CM (x), (23) 

being -f c the standard Dirac matrices, and ^> D v ^ = ^>D^ — (D^/)^f. Let us consider the FRW metric of the form 

ds 2 = dt 2 - a 2 (t) (^j + r 2 d9 2 + r 2 sin 2 (0)# 2 ) , (24) 
with p = 1 — «r 2 , and k = 0, 1, — 1. In the metric (24), the most natural choice for tetrads is 

e° = e^a; = dt, e 1 = e\dx Y = ^ydr, e 2 = e\dx 2 = a(t)rd9, e 3 = e\dx 3 = a(t)rsm(6)dr. (25) 
Using (19) and the definitions in (23) we have (see the Appendix) 

D = d + \[ li M^^ + \hulA—4 + \[^M t ^^^l * = 1,2,3; .7 = 2,3. (26) 

4 4 (XT' t: (XT' 

Since we choose diagonal tetrads every term eg with k = 1,2,3 is null. This implies that all the terms but the 
first arc null whatever the value of the commutators of the Dirac matrices. It is worth to stress that this result 
is independent of the choice of tetrads because the (0,0) component of the energy momentum tensor of a fluid is 
equivalent to the energy density only if the tetrads (or in general the chosen coordinates) are time-orthogonal as in 
(25). If we choose non time-orthogonal tetrads, i.e. eg not constant, T 00 does not represent the energy density because 
it acquires "pressure components" due to the different orientation of the tetrad. In our calculation these terms are 
the second, third and forth term of the (26). 

Thus the temporal component of the spinorial derivative in the FRW metric is just the standard time derivative: 

D = d . (27) 

This is not surprising if we consider the symmetries of the metric element (24). Thus, the (0,0) component of the 
stress energy tensor is 

Too - T^ lat . (28) 
This allows us to use T F lat to compute the cosmological constant. From Eq.(22) we thus obtain 

Too(x) = l - : (§„ l (x) 7 o < a > * m (x)) : (29) 

where : ... : denotes the customary normal ordering with respect to the mass vacuum in the flat space-time. 
In terms of the annihilation and creation operators of fields v\ and v 2 , T o is given by 



(30) 



and we note that Tqq is time independent. 

Next, our task is to compute the expectation value of Tqq in the flavor vacuum |0)^, which, as already recalled, 

is the one relevant to mixing and oscillations. Thus, the contribution (p™c) °f tnc neutrino mixing to the vacuum 
energy density is: 

f (o\j2 T w(mf = (p7:?)voo. (si) 

i 

We observe that within the above QFT formalism for neutrino mixing we have /(0|T O g |0) / = /(0(t)|T g |0(t)) / for 
any t. We then obtain 

f (0\J2 T oo(°)\°)f = J2 j ^ 3 k^(/(0|<X,,|0)/ + /(0|/3k! 4 /3LlO) / ). (32) 
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Since [1] 



f(OK\<Mf = f(0\Pl] t PlMf = sin 2 #|U k | 2 , (33) 

we get 

/<0|E T °o(°)l°>/ = 8sin ^ / ^Ki+wm) l^l 2 = (p™>oo, (34) 



i.e. 



(pZZ) = 327r 2 sin 2 / dk fc 2 ( W fc,i + w fe , 2 )|U k | 2 , (35) 
Jo 

where the cut-off if has been introduced. Eq. (35) is our result: it shows that the cosmological constant gets a non- 
zero contribution induced from the neutrino mixing. Notice that such a contribution is indeed zero in the no-mixing 
limit when the mixing angle 9 — and/or mi = m 2 (see Eq.(16)). It is to be remarked that the contribution is zero 
also in the limit of Vk — > 0, namely in the limit of the traditional phenomcnological mixing treatment. 

It is interesting to note that, for high momenta, the function |Vk| produces a drastic decrease in the degree 
of divergency of the above integral, in comparison with the case of a free field. Thus, if for example we chose 
K ^> sJm^rrH, we obtain: 

{pZc) oc sin 2 9 (m 2 - m x f K 2 , (36) 

whereas the usual zero-point energy contribution would be going like K 4 . 

Of course, we are not in a position to make our result independent on the cut-off choice. However, was not this 
our goal. What we have shown is that a nonzero contribution to the value of the cosmological constant may came 
from the mixing of the neutrinos. We have not solved the cosmological constant problem. Although it might be 
unsatisfactory from a general theoretical point of view, we may try to estimate the neutrino mixing contribution by 
making our choice for the cut-off. Since we are dealing with neutrino mixing, at a first trial it might be reasonable to 
chose the cut-off proportional to the natural scale we have in the mixing phenomenon, namely k 2 , ~ mim 2 [1]. 

With such a choice, using K ~ fc , mi = 7 x lCT 3 eU, m 2 = 5 x 10~ 2 eU, fc = 10~ 3 eV and sin 2 9 ~ 1 in Eq. (35), 
we obtain 

(p™T> = 1-3 x 10- 47 GeV 4 (37) 
Using Eq.(18), we have agrement with the upper bound given in Section I: 

A~10- 56 cm- 2 , (38) 
Another possible choice is the electroweak scale cut-off: K ~ lOOGeV. We then have 

(pZ7) = 1-5 x 10~ 15 GeV 4 (39) 

and 

A~l(r 24 cm- 2 , (40) 
which is, however, beyond the accepted upper bound. 



IV. CONCLUSIONS 



In this paper we have shown that the neutrino mixing can give rise to a nonzero contribution to the cosmological 
constant. We have shown that this contribution is of a different nature with respect to that given by the zero-point 
energy of free fields and we estimated it by using the cut-off given by the natural scale of the neutrino mixing 
phenomenon. The different origin of the mixing contribution also manifests in the different ultraviolet divergency 
order (quadratic rather than quartic, see Eq.(36)). The obtained value is consistent with the accepted upper bound 
for the value of A. On the contrary, by using as a cut-off the one related with the electroweak scale, the A value is 
greater than such upper bound. 
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It is worth to stress once more that the origin of the present contribution is completely different from that of 
the ordinary contribution to the vacuum zero energy of a massive spinor field. As we have shown, the effect we 
find is not originated from a radiative correction at some perturbative order [17]. Our effect is exact at any order. 
It comes from the property of QFT of being endowed with infinitely many representations of the canonical (anti- 
commutation relations in the infinite volume limit. Therefore, the new result we find is that it is the mixing 
phenomenon which provides such a vacuum energy contribution, and this is so since the field mixing involves unitary 
inequivalent representations. Indeed, as Eqs.(33)-(35) show the contribution vanishes as Vk — > 0, namely in the 
quantum mechanical limit where the representations of the (anti-) commutation relations are all each other unitarily 
equivalent. Our result thus discloses a new possible mechanism contributing to the cosmological constant value. 

As a final consideration, we observe that this effect could also be exploited in the issue of dark energy without 
introducing exotic fields like quintessence. In fact, neutrinos constitute a cosmic background of unclustered components 
whose mixing and oscillations could drive the observed accelerated expansion. 
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APPENDIX A 



For the sake of completeness, we derive the 0-th component spinorial derivative used to get the results (28). We 
start from the Cartan Equations (19) writing the tetrads one-forms in the FRW metric (24) as in (25). 
The exterior derivative of e a is 

de°=dtAdt = (Al) 

de 1 = - dt A dx 1 = - e° A e 1 (A2) 
p a 

de 2 = dr dt A d6 + a dr A d9 = - e° A e 2 + — e 1 A e 2 (A3) 

a ar 

de 3 = ar sin(6>) dt A d<p + a sin(0) dr A d<p + ar cos(#) d9 A d<p = 

= « e Ae 3 + ^ e l Ae 3 + ^) e2Ae 3. (A4) 

a ar ar 

Since the connection forms are antisymmetric in the tetradic (latin) indexes, Eq. (19) gives 

wJJ = (A5) 

w? = ~ e l = H e l i = 1,2,3 (A6) 

u) = — e? 3 = 2,3 (A7) 

J ar 

e 3 (A8) 

ar 

4 = 0. (A9) 



tan(0) 3 



Using the definition of spinorial derivative and spinorial connection given in (23), we have 

d„ = a„ + \hno] < + IbmA + \b2, 73] , (Aio) 

where we have used the antisymmetry of the commutators. The 0-th component of Eq.(AlO) gives Eq.(26). 



[1] M. Blasone and G. Vitiello, Annals Phys. 244 (1995) 283 [Erratum-ibid. 249 (1995) 363]. 



G 



[2] K. C. Hannabuss and D. C. Latimer, J. Phys. A 33 (2000) 1369; J. Phys. A 36 (2003) L69. 

[3] M. Binger and C. R. Ji, Phys. Rev. D 60 (1999) 056005; C. R. Ji and Y. Mishchenko, Phys. Rev. D 64 (2001) 076004; 

Phys. Rev. D 65 (2002) 096015. 
[4] K. Fujii, C. Habe and T. Yabuki, Phys. Rev. D 59 (1999) 113003 [Erratum-ibid. D 60 (1999) 099903]; Phys. Rev. D 64 

(2001) 013011; K. Fujii, C. Habe and M. Blasone, [hep-ph/02 12076]. 
[5] M. Blasone and G. Vitiello, Phys. Rev. D60 (1999) 111302; M. Blasone, P. Jizba and G. Vitiello, Phys. Lett. B 517 

(2001) 471; M. Blasone, P. P. Pacheco and H. W. Tseung, Phys. Rev. D 67 (2003) 073011. 
[6] M. Blasone, P. A. Henning and G. Vitiello, Phys. Lett. B 451 (1999) 140; 
[7] M. Blasone, A. Capolupo, O. Romei and G. Vitiello, Phys. Rev. D 63 (2001) 125015. 
[8] M. Blasone, A. Capolupo and G. Vitiello, Phys. Rev. D66 (2002) 025033. 

[9] N. Bahcall , J. P. Ostriker , S. Perlmutter , P. J. Steinhardt Science 284 (1999) 1481; S. Perlmutter et ai, Astrophys.J. 517 
(1999) 565. 

[10] V. Sahni, A. Starobinsky, Int.J.Mod.Phys. D9 (2000) 373. 

[11] Ya. B. Zeldovich, I. D. Novikov, Structure and evolution of the universe Moscow, Izdatel'stvo Nauka, (1975). 
[12] Ya.B. Zeldovich, Sov. Phys. Uspekhi 11, (1968) 381. 
[13] S. Weinberg, Rev. Mod. Phys. 61 (1989) 1. 

[14] M. Gell Mann, R.J. Oakes, B. Renner, Phys. Rev. 175 (1968) 2195; M.A. Shifman, A.I. Vainshtein, V.I. Zakharov, Nucl. 
Phys. B147 (1978) 385. 

[15] S. Capozziello, R. de Ritis, Phys. Lett. 203A, (1995) 283; Gen. Rel. & Grav. 29, (1997) 1425 ; S. Capozziello and R. de 

Ritis, Gen. Rel. Grav. 29 (1997) 1425. S. Capozziello, R. de Ritis, A. A. Marino, Gen. Rel. & Grav. 30, (1998) 1247. 
[16] S. M. Bilenky and B. Pontecorvo, Phys. Rep. 41 (1978) 225. 
[17] S. Coleman, E. Weinberg, Phys. Rev. D 7 (1973) 6. 



7 



